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A sensitivity of the VEPP-2000 e + e _ collider in a search for the rare decay r/ —> e + e~ has been studied. 
The inverse reaction e + e“ -A r] is proposed for this search. We have analyzed a data sample with an integrated 
luminosity of 108 nb -1 collected with the SND detector in the center-of-mass energy range 520 — 580 MeV 
and found no background events for the reaction e + e“ -A r/ in the decay mode rj — > 7 r 0 7r°7r°. In the absence of 
background, a sensitivity to Bp —> e + e“) of 1CP 6 can be reached during two weeks of VEPP-2000 operation. 
Such a sensitivity is better than the current upper limit on Bp —¥ e + e _ ) by a factor of 2.3. 
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Figure 1. The Feynman diagram for ri —> e + e . 

Decays of pseudoscalar mesons to the pair of lep¬ 
tons P -A l + l~ are rare. In the Standard Model (SM) 
these decays proceed through the two-photon interme¬ 
diate state as shown in Fig. |T] and therefore are sup¬ 
pressed by a factor of ~ a 2 comparing with the P —> 77 
decays, where a is the fine structure constant. An ad¬ 
ditional suppression of ( mi/mp ) 2 arises from the ap¬ 
proximate hclicity conservation, where to; and mp are 
the lepton and meson masses, respectively. So, the 
P -A l + l~ partial width is less than the two-photon 
width r(P —> 77) by a factor of ~ a 2 (mi/nip) 2 . The 
low probability makes these decays sensitive to possible 
contributions of new physics beyond the SM Q] 2], In 
the SM, the knowledge of the transition meson-photon 
form factor F(q\, q %) for the 7*7* —1 P vertex is needed 
for a r(P -A l + l~) calculation, where qf and q 2 are four- 
momenta squared of virtual photons in Fig. [L] The real 
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and imaginary parts of the decay amplitude are usu¬ 
ally calculated separately. The imaginary part 3(A) 2 
is proportional to the form-factor value |F(0,0)| 2 at 
qf = q 2 = 0 and, consequently, can be calculated us¬ 
ing the known width of the P —> 77 decay. Since 
|A| 2 > 3(A) 2 , a model-independent lower boundary 
(unitary bound) on the P —1 l + l~ width can be ob¬ 
tained from r(P — > 77) [3]. The real part of ampli¬ 
tude cannot be calculated in a mo del-independent way. 
In Ref. [4] it is shown that the real part depends on 
the integral over q 2 < 0 of the form factor in symmet¬ 
ric kinematics, F(q 2 ,q 2 ). The ranges of predictions for 
the P -a l + l~ branching fractions obtained in different 
form-factors models 0E| are listed in Table I. For com¬ 
parison, the last column of Table I contains the current 
experimental values of the branching fractions. Cur¬ 
rently, only two of the five decays are measured and for 
two more upper limits on the width exist. The value 
of Bp -A pPp~) agrees with the prediction but the 
measurement accuracy is low. The more precise value 
of Bp 0 -a e + e ~) differs from the prediction by about 
three standard deviations. 

It is clear that the current experimental situation 
requires improvements in accuracy for the n° -A e + e~ 
and 77 -A gPp~ decays as well as measurements of the 
other three decays. Such measurements are planned at 
the BES-III H2], KLOE-2 QA], Crystal Ball QH and 
WASA [T3, TBj detectors in the near future. 

The upper limit on Bp' -A e + e~) has been recently 
set in experiments with the CMD-3 [TT] and SND [TO] 
at VEPP-2000, where the inverse reaction e + e“ -A rj 1 
was used. In this paper we propose the same method 
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Table I. The unitary bounds, theoretical predictions in units of the unitary bound E1IS, and experimental values for 
the P —> l + l~ branching fractions. 


B(p->i+r) 

Unitary bound (U) 

Theory 

Experiment 

B( 7 T° -s- e+e") x 10 8 

4.69 

6.23-6.38 

7.49 ±0.38 0 

B(rj —*■ e+e - ) x 10 9 

1.78 

4.60 - 5.24 

< 2300 [5] 

B(rj n + n~) x 10 7 

4.36 

4.64-5.12 

5.8 ±0.8 0 

B(r)' —> e+e - ) x 10 10 

0.36 

1.15-1.86 

< 56 EDI EH 

B(rj' x 10 7 

1.35 

1.14-1.36 

- 


for a search for the r) —> e+e - decay and perform its 
feasibility study. 

Data used in this paper were collected with the SND 
detector at the VEPP-2000 [17] collider in 2013. VEPP- 
2000 is designed for a study of e+e - annihilation at the 
center-of-mass (c.m.) energy (E) from 160 MeV up to 2 
GeV. There are two detectors at the collider, SND and 
CMD-3, which collect data simultaneously. At the mo¬ 
ment of data taking the VEPP-2000 accelerator complex 
consisted of the 3 MeV electron linac ILU, the 250 MeV 
pulsed synchrotron B-3M, the 900 MeV booster storage 
ring BEP, and the collider storage ring VEPP-2000. 

Currently the complex VEPP-2000 is being up¬ 
graded. The maximum energy of BEP will be in¬ 
creased up to 1 GeV, and the injection system will be 
changed. Electrons and positrons will be transported to 
BEP from the VEPP-5 injection complex |18| through 
a 250 m beamline. The upgrade allows to increase the 
VEPP-2000 luminosity at maximum energy up to 10 32 
cm- 2 sec _1 and should result in a more stable operation 
of the accelerator complex. 

For the current study, VEPP-2000 parameters at 
c.m. energy close to m^c 2 = 548.862 ± 0.018 MeV [9] 
such as luminosity, accuracy of the energy setting, en¬ 
ergy spread, are important. In 2013 SND did not record 
data exactly at this energy. Therefore, we analyze data 
from four energy points near m v c 2 , with c.m. energies of 
520, 540, 560, and 580 MeV. The integrated luminosity 
collected at these energy points is measured using the 
reaction e+e - —> 77 to be 108.1±2.0 nb _1 . The average 
luminosity during data taking varied from 0.26 x 10 3 ° 
cm _ 2 sec _1 at 520 MeV to 0.73 x 10 3 ° cm _ 2 sec _1 at 580 
MeV. Interpolating the energy dependence we estimate 
that the average luminosity expected at E = m v c 2 is 
0.34 x 10 3 ° cm _ 2 sec _1 . 

The width of the ?? resonance, = 1.31 ± 0.05 
keV [5], is much less than the collider c.m. energy 
spread <rg. This means that the value of erg is cru¬ 
cial for the proposed search for the e+e - —> r/ reaction. 
At E = m v c 2 the energy spread is determined by the 


Touschek effect. It is uniquely related to the RMS of the 
longitudinal distribution of the e+e - interaction point 
<Jz EH: 

erg = 4.05a z ^/v r cav £ ; bSin[arccos(63.2E;^/V cav )], (1) 

where erg is measured in keV, erg in mm, the RF cavity 
voltage 14av in kV, and the beam energy E^ = E/2 in 
GeV. The value of erg is measured using detected events 
of the elastic scattering e+e - —> e+e - . At E = 540 
MeV data taking was carried out with several values of 
the RF cavity voltage ranged from 18 kV to 36 kV. The 
values of erg varied from 16 to 12 mm, while the values 
of erg from 120 to 160 keV. Below we will use the value 
of erg = 150 keV as an estimate of the energy spread at 
E = m v c 2 . 

In the proposed experiment the collider energy must 
be set and monitored with an accuracy better than erg. 
This is provided by the beam-energy-measurement sys¬ 
tem based on measuring energies of Compton back- 
scattered laser photons [15] , The accuracy of this 
method is about 60 keV. The system monitors beam 
energy continuously during data taking; the single mea¬ 
surement duration at Em 550 MeV is about one hour. 

The search for e+e“ —>7 events will be performed 
with the SND detector [201 1211 [22] [23]. It is a non¬ 
magnetic detector. The main part of the detector is 
a three-layer spherical electromagnetic calorimeter con¬ 
sisting of 1640 Nal(Tl) crystals. The calorimeter covers 
a solid angle of about 90% of 47 t. The energy and angu¬ 
lar resolutions for photons with energy U 7 are described 
by the following formulae 


ag 7 /U 7 = 4.2%/^/u 7 (GeV), 

( 2 ) 

&e,(t> = 0.82°/i/ .E(GeV). 

(3) 


Directions of charged particles are measured in a nine- 
layer drift chamber. The calorimeter is surrounded by 
an iron absorber and a muon detector. In the proposed 
search the muon detector is used to veto cosmic rays. 
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The energy dependence of the Born cross section for 
the reaction e + e _ —> 77 is described by the Breit-Wigner 
formula: 


CTo = 


e + e 




(ml — E 2 ) 2 + mlVl 


(4) 


In experiment, it is necessary to take into account ra¬ 
diative corrections arising, for example, from radiation 
of extra photons from the initial state. To do this, the 
cross section is convolved with the so-called radiator 
function W{s,x) [23j 1251 


o(s) = f W(x, s)cto(s(1 — x))dx, (5) 

J o 


where s = E 2 . The upper limit of integration in Eq.® 
depends on the decay mode and equals unity for the rj 
decay into two photons and 1 — (3m w o) 2 /s for the de¬ 
cay into 37 t°. The theoretical accuracy of the corrected 
cross section ® is better than 1% [23J 25). For the 
unitary bound B(r] —> e + e _ ) = 1.78 x 10 “ 9 the Born 
cross section in the resonance maximum is oq = 29 pb. 
The radiative corrections lead to decrease of this cross 
section down to o = 14 pb. 

To take into account the collider energy spread, the 
cross section ® should be convolved with the Gaussian 
function describing the distribution of the integrated lu¬ 
minosity over energy 


+oo 


-'exp 


(E 0 ) = 


\Z2ttoe 


(E-EnT 

e 2ct e <j(E)dE, 


( 6 ) 


where Eq is the average beam energy. For oe = 150 
keV, E 0 = m v c 2 , and B(r/ e + e~) = 1.78 x 10 -9 the 
visible cross section is 


ofx P (wr,c 2 ) = 105 ± 11 fb. (7) 

Thus, the energy spread allows to use only cr e xp/o’ ~ 
1/140 of collected integrated luminosity. The uncer¬ 
tainty on <jg/ p (m^c 2 ) is due to the uncertainty of the 
beam-energy measurement (60 keV). 

The main 77 decay modes are r) —> 77 (39.4%), 
77 —» 7 r° 7 r° 7 r 0 (32.7%) and 77 —> 7 r + 7 r _ 7 r° (22.9 %). The 
numbers in parentheses are the branching fractions [9]. 
In the 77 —> 77 mode the process e“*"e - —> 77 cannot be 
separated from the QED process e + e _ —> 77 . 

The most suitable 77 decay mode for the search for 
the e + e _ —> 77 reaction at SND is 77 —> n a n a Tr° —> 67 , for 
which physical background is small. The main source 
of background is cosmic-ray events. For the search for 
e + e _ —> 77 , events with six or more detected photons 
and with the energy deposition in the calorimeter larger 



Figure 2. The distribution of y 2 of the kinematic fit 
for simulated e + e” —>• 77 events passing the selection 
criteria for six-photon events. 


than 0.6 E are selected. Background from events with 
charged particles is rejected by the selection condition 
that the number of fired wires in the drift chamber is 
less than four. Cosmic-ray background is suppressed by 
the veto from the muon detector. 

For events passing the preliminary selection, a kine¬ 
matic fit to the e + e _ —> 7r 0 7r°7r° —> 67 hypothesis is 
performed with a requirement of total energy and mo¬ 
mentum conservations and a condition that invariant 
masses of three photon pairs are equal to 7 r° mass. The 
quality of the kinematic fit is characterized by the pa¬ 
rameter y 2 . During the fit all possible combinations of 
photons are checked and a combination with the small¬ 
est x 2 value is selected. The x 2 distribution for sim¬ 
ulated events of the process e + e _ —> 77 passed the se¬ 
lection criteria for six-photon events is shown in Fig [2] 
The condition x 2 < 100 is used to select 77 candidates. 

The detection efficiency for e + e _ —>77 events is ob¬ 
tained using simulation to be £ = (12.5 ± 0.6)%. The 
quoted error is systematic. It is estimated using results 
of Ref. [2Tj, where comparison of data and simulated x 2 
distributions was performed for five-photon events from 
the process e + e” —>• ujtt 0 — > 7r°7r°7. 

No events satisfying the six-photon selection criteria 
are found in the data sample with an integrated lumi¬ 
nosity of 108.1 ± 2.0 nb -1 recorded in the energy points 
with E = 520, 540, 560 and 580 MeV. 

Main sources of physical background for the decay 
mode 77 —> TT°TT°ir 0 are five-photon events of the pro¬ 
cesses e + e _ —> 7r°7r°7 and e + e“ —> 5y with a spurious 
photon from beam-generated background or splitting of 
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electromagnetic showers in the calorimeter. The prob¬ 
ability to find a spurious photon in an event may reach 
up to 5%. This number is used below for background 
estimation. The cross section for e + e _ —> 7 r° 7 r °7 at 
E w 550 MeV is determined mainly by the transition 
e + e _ —> p( 770) —> 7 r° 7 r °7 and is estimated to be 2 ± 1 
pb. The quoted uncertainty is due to an unknown mech¬ 
anism of the the p —> 7 r° 7 r °7 decay |28| . The detection 
efficiency for e + e _ — > 7 r° 7 r °7 events with the six-photon 
selection criteria is about 1%. The cross section for 
e + e _ —> 5y with detection of all five photons is cal¬ 
culated using the CompHep program [29] to be about 
4 pb. The requirement of an extra photon and con¬ 
dition x 2 < 100 decrease the cross section by a fac¬ 
tor of about 1/200. Thus, the total cross section of 
background processes with the selection criteria for the 
decay mode 77 —> 7 r 0 7 r° 7 r° is about 0.04 pb and negligi¬ 
ble at integrated luminosities, which can be collected at 
VEPP-2000. 

The second decay mode suitable to search for 
e + e _ — > 77 events is 77 —>• 7 r + 7 r _ 7 r°. For this mode 
there is a background from the nonresonant process 
e + e _ — > tt + tt~tt°. There are no data on the e + e _ —> 
cross section at Em 550 MeV. To estimate it, 
we use the results of Ref. [30], where the cross section 
e + e _ was measured by the SND detector 

at the VEPP-2M collider in the energy region of the 
w(782) and p(770) resonances. In Ref. [30] the data 
were fitted in several models, which at E = 550 MeV 
give the cross section value ranged from 6 to 16 pb. Such 
cross section values correspond to branching fractions 
B(i 7 — > e + e~) = (0.4 — 1.2) x 10 -6 . Thus, at the level 
of sensitivity to B{rj —>■ e + e~ ) of 10 -6 the measurement 
of the e + e _ — > out of the 77 -meson resonance 

may be required. Other background sources in the mode 
77 TT + Tr~n° are the QED processes e + e _ —> e + e~j( , y) 
and e + e _ — > p + In a data sample collected with 
the SND detector at E = 520, 540, 560, and 580 MeV, 
we cannot suppress background from these processes to 
the level reached in the mode e + e“ —> 7r 0 7r°7r°. 

The measured cross section e“*"e _ —> 77 is determined 
as 

N s 

CTexp = eL’ ( 8 ) 

where N s is the number of selected e + e _ —> 77 events, e 
is the detection efficiency, and L is the integrated lumi¬ 
nosity. 

From our result that no background events for the 
decay mode 77 —> 7 r° 7 r 0 7 r° are found in the data sample 
with an integrated luminosity of 108.1 nb -1 we estimate 
the sensitivity of the e + e _ —> 77 cross section measure¬ 


ment. The upper limit on the cross section at the 90% 
confidence level (CL) [31; corresponding to N s = 0 is 


2.3 

CTexp < 0.125- 108 


170 pb. 


(9) 


Comparing this limit with the cross section 0 calcu¬ 
lated for the unitary bound B (77 —> e + e _ ) = 1.78 x 10 -9 , 
we estimate the sensitivity to the search for the decay 
77 —> e + e~ with an integrated luminosity of 108 nb -1 to 
be 

£(77 e+e“) < 2.9 x 10 " 6 (10) 

at 90% CL. This result is close to the upper limit 
B(i 7 —t e + e _ ) < 2.3 x 10 -6 set recently in the HADES 
experiment 0 , in which 77 mesons were produced in 
proton-nucleon collisions. With a VEPP-2000 luminos¬ 
ity of 0.34 x 10 3 O cm _ 2 sec _1 the current upper limit can 
be reached in a week of data taking. In two weeks a 
sensitivity at the level of 10~ 6 can be reached. 

In conclusion, we have analyzed the data sample 
with an integrated luminosity of 108 nb -1 collected with 
the SND detector at the VEPP-2000 e + e _ collider at 
c.m. energies 520, 540, 560, and 580 MeV and found 
no background for the process e + e _ —> 77 in the de¬ 
cay mode 77 —> 7r 0 7r°7r°. In the absence of background, 
data with an integrated luminosity of 324 nb -1 provide 
a sensitivity of 10 ~ 6 for U (?7 —> e + e~). Such data may 
be accumulated in two weeks of VEPP-2000 operation. 
The sensitivity of 10 ~ 6 is 2.3 times better than the cur¬ 
rent upper limit on B(r] —> e + e - ). 

Data at E = m^c 2 will be recorded during the 
experiment on a measurement of hadronic cross sec¬ 
tions below 1 GeV planned at VEPP-2000. Two detec¬ 
tors, SND and CMD-3, will collect data simultaneously. 
Their results may be combined to improve sensitivity to 
B(rj —> e + e _ ) by a factor of about 2. 
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